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ABSTRACT 

We argue that gravitational instability of typical protostellar disks is not a viable mechanism 
for the fragmentation into multiple systems - binary stars, brown dwarf companions, or gas 
giant planets - except at periods above roughly 20,000 years. Our conclusion is based on a 
comparison between prior numerical work on disk self-gravity by Gammie (2001) with our own 
analytical models for the dynamical and thermal state of protostellar disks. For this purpose 
we first develop a simple theory for the initial conditions of low-mass star formation, accounting 
for the effect of turbulence on the characteristic mass, accretion rate, and angular momentum 
of collaping cores. We also present formulae for the probability distribution of these quantities 
for the case of homogenous Gaussian turbulence. However, our conclusions are not sensitive to 
this parameterization. 

Second, we examine the criterion for fragmentation to occur during star formation, con- 
centrating on the self-gravitational instabilities of protostellar accretion disks in their main 
accretion phase. Self-gravitational instabilities are strongly dependent on the thermal state of 
the disk, and we find that the combination of viscous heating and stellar irradiation quenches 
fragmentation due to Toomre's local instability. Simulations by Matsumoto & Hanawa (2003), 
which do not include detailed thermal evolution, predict fragmentation in an early phase of 
collapse. But, fragments born in this phase are on tight orbits and are likely to merge later 
due to disk accretion. Global instability of the disk may be required to process mass supply, 
but this is also unlikely to produce fragments. We conclude that numerical simulations which 
predict brown dwarf formation by disk fragmentation, but which do not account for irradiation, 
are unrealistic. Our findings help to explain the dearth of substellar companions to stellar type 
stars: the brown dwarf desert. 



1. INTRODUCTION 

At birth, stars accumulate their material through 
disks that are well known to be dense and mas- 
sive. Both of these properties make protostel- 
lar accretion dis ks susceptible to self-gravitational 
instability (e.g., iToomrd . 



19641 ). Indeed, simula- 



tions of pr otostellar d isk accretion and evolution 
like those of iLin &: Pri ngle (1990) typically find that 
they approach or cross the threshold for instabil- 
ity, either because of a high disk density or be- 
cause of a finite disk mass. When this occurs, 
disk motions caused by self-gravity become impor- 
tant and perhaps d ominan t sources of angular mo- 
ment u ni transport (iPaczvnski 1978 ; Lin &: Pringle 
19871 : iPapaloizou k. SavonilT Il99ll: lAdams et aO 



1989*; 'Tohhue* Il994l : iLaughlin fc Bodenheimeil Il994l : 
Gammie 2001). Indeed, they may be required for 
disks to process mass accretion. 

Disk self-gravity is of interest not only as a trans- 



port mechanism, but also because it may cause the 
disk to fragment into self-gravitating bodies that can 
become stellar, substellar, or possibly planetary com- 
panions to the star being formed. 

If fragments form sufficiently early and survive as 
most of the stellar mass is accumulated, they stand 
a chance of acq uiring mass comparable to that of the 
prim ary object (lBonnelll994l :l IMatsumoto &: Hanawal 
200^. Alternatively, fragments that survive but fail 
to accumulate material may wind up with brown 
dwarf or planetary masses at the end of accre- 
tion, depending partly on where in the disk they 
form. The possibility that gas giant planets form 
fro m fragmen t ation of a gaseous disk was suggested 
by 'Cam eronI ( 19781 ) and advocated by Boss (e.g., 
Bosail997^. I t has been investigated numerically by 
Maver et al. ( 20021 ). for example. However, numer- 



ical investigations are difficult, both because of the 
high resolution required to attain convergenc e in self 



gravitating simulations of cooling gas ((Truelove et al. 
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^93), and because of the imp ortance of thermal 
evolution to d isk fragmenta tion (iNelson et al.ll200Cll : 
GammidEnnil : iPickett et al1l2nn,'j : lE,afikovll2nn4 ) . We 

return to disks' thermal evolution below. 

The possibility that disk instability can produce 
stellar or substellar companions makes it relevant to 
the problem of binary and multiple star formation, 
and to the origin of the stellar initial mass func- 
tion. The binary problem is reasonably well-posed, 
since the p restellar cores have been mapped (e.g., 
Benson fc Mvers 1989. ) and binary stars are also well 
studied. Inde ed, binaries' angular momenta are sim- 
ilar t o cores' ( Duauennov Mavojll99il : IBoss et al.l 
2000l ). a f act that exp lains some trends in binary 
properties ( Fished 2004 ) . However the mechanism by 
which a core fragments is currently unclear. It is not 
understood, for instance, why some stars - albeit less 
t han half - form on their own. 

iTohlind ((20o2) has reviewed theories for the origin 
of binary stars, and identified three leading mech- 
anisms: "prompt" fragmentation in a sheet formed 
from the quasi-homologous collapse of a cold, rotat- 
ing cloud; fission of a rapidly rotating and contracting 
protostar, a nd self -gravitational instabilities of accre- 
tion disks. ( Hovld 's 1953 proposal of opacity limited 
fragmentation during collapse appears unlikely, due 
to pressure gradients caused by central concentration 
of the initial state.) One should add to Tohline's list 
turbulent fragmentation, in which nonlinear perturba- 
tio ns in the ini tial state engender fragments directly; 
see IPadoan fc N ordlund (2002), Klein et al. (2003), 
KlessenI l^lOol ). and lBate et al.l (I2003D . 

This paper will consider axisymmetric models of 
low-mass star formation. Given the turbulent charac- 
ter of the molecular clouds from which stars form, one 
might wonder how relevant axisymmetric models can 
possibly be to stellar fragmentation. However, the en- 
ergy in this t urbulence is concentrated in its longest 
wavelengths The direct progenitors of 

single and binary stars are the dense molecular cores 
which are supported significantly by thermal pres- 
sure, unlike the larger structures that contain them - 
their parent cloud or the cloud substructures known 
as clumps - which are sup ported entirely bv t urbu- 
lence and mag netic fields (|McKee et aD " ll993r ). On 
the scale of an individual prestellar core, perturba- 
tions with longer wavelength appear as translational 
motion, shear, and overall compression. All of these 
can be treated with an axisymmetric model in the 
center of mass frame. Motions with smaller wave- 
lengths constitute turbulence within the core and 
contribute to turbulent fragmentation but are weaker 
in general and are suppresse d by ion-neutral friction 
( Zweibel fc .TosafatssonI 198^ . For these reasons, and 
for analytical convenience, we will concentrate on ax- 



isymmetric initial conditions. 

In this paper we develop a simple prescription for 
the initial conditions of star formation and examine 
the susceptibility of the predicted protostellar accre- 
tion disks to several modes of gravitational fragmen- 
tation. We begin by modeling the effect of molecular 
cloud turbulence on the masses, radii, and rotation 
rates of the prestellar molecular cores (§121). We then 
address the "prompt" instability that occurs early in 
the collapse (§151) , using criter ia develope d bv t he nu- 
merical survey of lMatsumoto fc Hanawa (f2003l ) . This 
analysis suggests that fragmentation occurs quite fre- 
quently in the early collapse. However, we argue that 
this may be an artifact of approximations made in 
the thermal evolution of core material, and that the 
fragments so formed are destined to merge during 
subsequent accretion. 

We t h en tu rn (§ ^ to fragmentation due to 
iToomrel (|l964l Vs instability. Since this is governed 
by the local density, we will call it "local fragmen- 
tation". Here we must discriminate between condi- 
tions in which Toomre's instability saturates, i nduc- 
ing angular n i qment u m transport dP aczvnski' '1978'; 
T;in Pringlel Il987l: IPanaloizoii Savoniie 1991: 



Laughlin &: Bodenheimeill994l : lLaughlin fc Rozvczkal 



1990), and conditions that l ead to lo c al fra gmen- 
tation. For this we rely on lOammid (|200lh . who 
identifies the boundary between saturated "gravi- 
toturbulence" and fragmentation in his simulations 
of razor-thin accretion disks. Although protostellar 
disks are hardly razor-thin, Gammie's criterion is a 
useful litmus test for fragmentation. It demonstrates 
that the outcome of gravitational instability is sen- 
sitive to the thermal evolution of the disk midplane, 
which we model first by treating only viscous heating 
alone in ^4.1.11 This analysis indicates the possibility 
of fragments forming with periods greater than about 
460 years. However when we account for irradiation 
from the central star in ^4.2| we will demonstrate 
that this quenches local fragmentation in the main 
accretion phase, at least for periods up to ~ 1.2 x 10^ 
years. 

Lastly we treat instabi litv t hat arises from the fi - 
nite disk mass (Adams et al. 1989: Tohlinej ^13) • 
This sets in if othe r transport m echanisms fail to 
fiush the disk, and iBonnelll ( 19941 ^ suggests it it may 
lead to binary fragmentation. Howe ver, current ev 



idence from numerical simulations ( Laughlin et al 



Il998l : iLaughlin fc Rozvcz ka"l99eh shows that the in 
stability would always redistribute the disk material 
to quench itself, and the fragmentation would never 
occur via this channel. 

We draw conclusions in ^6.31 about the implica- 
tions of this work for understanding the brown dwarf 
desert, i.e., the observed underrepresentation of sub- 
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solar companions to solar type stars. 

2. STAR FORMATION: INITIAL CONDITIONS 

Prior to their collapse, molecular cores are partially 
supported by thermal pressure and are confined by 
the hydrostatic pressure of their parent cloud and 
permeated by its turbulent motions. Combined with 
the scalings known to apply to molecular cloud tur- 
bulence, these statements suffice to predict, if only 
cru dely, th e initial conditions for star formation. See 
also lMcKee (1999) and McKee k Tan (2003). 

First, consider the structure of cores prior to their 
collapse. In the absence of non-thermal pressure from 
magnetic fields, tu rbulence , or r otation, these are 
Bonnor-Ebert (BE; iBonnod 119561 ) spheres, confined 
by some external pressure P. They exhibit a sequence 
of increasing self-gravity and central concentration, 
leading up to a critical state whose mass and radius 
satisfy 



Mbe = 1.18 



a. 



th 



G3/2pl/2 



and 



Mbe 
-Rbe 



2.43-^, (1) 



where fith isothermal sound speed. Molecular line 
cooling and dust radiation typically conspire to main- 
tain (Tth — 0.2 km/s, correspo nding to temperature s 
T ~ 10 K for molecular gas (iFuller fc Mverslll993^ : 
however one can imagine other possibilities, such as 
star formation in strongly heated or metal-free gas, 
in which T could be significantly higher. 

Because non-thermal forms of pressure also con- 
tribute, it is reasonable to replace fith with an ef- 
fective sound speed Ces derived below. The Bonnor- 
Ebert parameters Mbe and i?BE provide fiducial val- 
ues for the core mass and radius, although the actual 
values (Mc and Rc, say) differ because CeS > fth- We 
show below that thermal and nonthermal contribu- 
tions are compara.ble in thes e objects. 

ICo odm an et al. ( ^9^) and Burkert &: Bodenheimei] 
(|2000|) have shown that the observed rotation rates 
of molecular cores are consistent with the turbulent 
velocity fields of their parent clouds, if extrapolated 
to the sizes of individual cores. It follows that ctnt 
for an individual core can be predicted the same way. 

Unlike cores, molecular clouds are supported by a 
triumvirate of turbulent velocity, turbulent magnetic 
fields (which together compose unt) and magnetic 
field pressure, with ath being entirely negligible on 
scales above Rc- Clouds are famously repo rted to 
obey a line width-siz e relation, noted first by iLarsoiil 
(|l98ll ) and refined bv lSolomon et al.l ((l2l3), in which 



cTNTiRf ^ 0.70GSciiZ, 



(2) 



corre sponding to a state of strong self-gravity 

199i). In eq. @ , R refers to the scale on which ctnt is 

measured - on the core scale, cjnt would be c7nt(-Rc)- 



The column densi ty of the paren t cloud or molec- 



ular clump is Sci; ISolomon et al. (|l987t ) find 



us- 
ing a virial analysis, that for giant molecular clouds 
is remarkably constant around the value 165 
Mq /p c^ . However , Sni does depend o n cori text : in 
M33, for example, lEosolowskv et af] (l2003l) derive 
Sci ~ 120M^ /nc^: in the LMC. IPak et aJ.1 ()l99Sl 'l 
derive Ed ~ TOOMq/pc^ (ifc0.3d ex): and in Galac- 



tic clumps forming massive stars ( Plume et al. 19971 : 



McKee T^ '2003'^ or in ext ragalacti c starburst re 
gions (e.g.. ^Sarge nt fc Scoville 1991i : IScoville et al 
1991), Sci ~ 3000 - 5OOOM0/pc2. 



The confining pressure for a core must derive from 
the mean hyd rostatic pressu re of the parent clump or 
cloud; this is (|McKeelll999l 'l 



P 



O.Sl/pGS^i. 



(3) 



where we have included fp to account for variations 
in clump or cloud structure (0.51 refers to a spherical 
cloud with p (X 1/r and no embedded stars) and for 
fluctuations around the mean value due to turbulence 
or location within the cloud. With equations this 
implies 



Mbe = 1.65 



'th 



Rbe = 0.68 



'th 



0.70 165Mopc-^ ^2 

-jT/2 ^flf.ioM©, 

0.034 165M0PC-2 



1/2 



Tcs,wpc (4) 



Ceff) 



where Tefj^io is the temperature for which ath 
measured in units of 10 K. 

We wish to highlight two points about Bonnor- 
Ebert spheres confined by the ambient pressure 
(equations P and [3]) and bathed in turbulence de- 
scribed by equation First, note that the column 
density of the core is similar to that of the parent 
cloud or clump: 



Mbe 



1.1/V2S 



cl, 



(5) 



independently of ath- This results from the fact that 
both the core and the region in which it is embed- 
ded are self-gravitating objects whose pressures follow 
equation ©, and also from the fact that the core's 
mean pressure is prop ortional to its surface pressure. 
McKee fc T^ (|200i) find a relation quite similar to 
eq. (0) for turbulence-supported cores. 
Second, compare ctnt(-Rbe) to ath- 



(Tnt{Rbe) 0.69 



c^th 



fi 



1/4- 



(6) 
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independently of Sd. Equation © implies that any 
cloud substructure whose column density is compa- 
rable to Sci (i.e, any self-gravitating substructure) is 
suffused with turbulence at the virial level. Extrap- 
olated to the scale of a Bonnor-Ebert sphere, which 
is assumed to be supported by thermal pressure, this 
implies ctnt — cth as in equation ®. Heuristically, 
the hydrostatic cloud or clump pressure P, which is 
nonthermal on scales above Rbe, must be continuous 
with the core pressure, which is largely thermal. 

Equation ©, in turn, implies that nonthermal mo- 
tions must be included in the support of cores. To be 
specific, we adopt ^ 



Ceff 



0-t\ + 1.06(T^T 



(7) 

motivated by McKee &: HollimanI (|l999l Vs analy- 
sis of the mass of regions supported by Alfven 
waves (treated in t he WKB approximation: 
Mr.Kee Ri. Zweihell Il995h . With this, equation ® 



is solved if 

0.98 

<TNT — 



1/2 



2.02 



ft 



1/2 



T, 



(cores) (8) 



which are power-law fits to the algebraic solution 
(good to 20% for 0.5 < fp < 3). Cores are there- 
fore somewhat larger than objects lacking turbulent 
support: 



Mr 



Rr 



4-1 

-^Mbe 
Jp 



2.0 

Jp 



Rbe 



2.8 165M0PC- 
0.069 165M0PC-2 



Tio pc (9) 



These values are more consistent with observations 
than i?BE and Mbej which points to the presence 
of turbulent support. The fact t hat Mr exceeds the 



in comparison to the cloud or clump mean pressure. 
This is a product of several factors: the distribution 
of pressure within a cloud; the distribution of pres- 
sure with time due to turbulent fluctuations; and the 
likelihood of forming an unstable core, which presum- 
ably biases the distribution to higher pressures. On 
the other hand, the fp distribution is strongly con- 
strained if one posits that stars always form from 
marginally stable, thermally-supported cores whose 
masses are given by equation Q. In that case, the 
stellar initial mass function is a direct consequence 

—3/2 

of the core mass distribution: Mc oc efp . Con- 
versely, the distribution of fp is empirically deter- 
mined, in such a theory, by the distribution of s tel- 
lar (a nd binary) masses. For instance, Padoan e t al.l 
(1993) attribute the stellar initial mass function to 
the statistical distribution of pressure in simulations 
of supersonic, isothermal, unmagnetized turbulence. 
We will not adopt this strategy for the distribution of 
/p, however, for two reasons. First, such a theory pre- 
dicts that the most massive stars form in the lowest- 
density regions, which contradicts observations. Sec- 
ond, cores that create massive stars are significantly 
more mas sive than Mrf, and therefore highly turbu- 
lent (e.g.. lMcKee TanlEofii ^. We reiterate that we 
consider only the typical, low-mass cores, for which 
thermal support is naturally significant. 

We have not considered static magnetic fields in 
our treatment of cores; we comment on this in ^2.41 



Core Rotational Properties 
Rurkert Rodenheimer] (|2nOft l and 



mean stellar mass is no concern: Matzner fc: McKe^ 
({2000') have shown that low-mass stars form at an 
efficiency e ~ 25% — 80%, due to the action of proto- 
stellar winds. Bear in mind, moreover, that there is 
a mass spectrum of cores going to larger and smaller 
masses, which generally res embles in its shap e the 
stellar initial mass function ( Motte et al.lll99g ). 

It is also worth stressing that the objects we call 
"cores" are only the smallest and most thermally- 
supported self-gravitating regions identifiable within 
molecular clouds. The clumps that form clusters of 
stars, and clouds themselves, have ctnt ^ Cth and 
masses much greater than Mc- 

We have not yet specified the distribution of the 
variable fp that determines cores' bounding pressures 

(|7|l depends on how Ccff is used 



2.1. 

We follow 

Fished (jsQOJ in assuming that the velocity field is 
homogeneous and Gaussian random with a spectrum 
consistent with equation These authors also as- 
sume the components of v to be uncorrelated, reduc- 
ing the mean vorticity relative to an incompressible 
velocity field of the same a. We shall first adopt 
the assumption of uncorrelated velocity components, 
then consider how the result should be modified to 
account for gas pressure. 

Recause both ctnt and the specific angular momen- 
tum j derive from v, one expects j oc i?cO'NT(-Rc)) 
with statistical fluctuations about this scaling. This 
expectation is borne out in Appendix where we 
compute the full distribution of j. Equations ()A19|) 
and ()A20|) imply, for a velocity field with uncorrelated 
components. 



J 



0.34/,ii,aNT 



2.3 X 



10' V, 



U65Mqpc- 



1/2 



Rr 



0.1 pc 



3/2 



cm 



^The coefficient of ct^t in ^l- 



For Alfven waves treated as a gas with adiabatic index 3/2 
and polytropic index 1/2, this coefficient is 3/2 for P = pcla JMcKee fc ZweibelllTOQl) . 2/3 for A/be = l.lScl^ /[G^''^ p{eAgef^] 
IpcKee fc HollimaiL.199a) . and 1.06 for Mbe = 1.18Ceff/(G^''^P^^^). 



5 



1.3 X 10 



21 



fj ^3/2 165Mqpc- 



3/2^10 



(10) 



where the last hne uses Rc from equation Q. 

Since we have assumed v to be scale free, fj is 
chosen from a distribution that depends only on the 
internal structure of the core, the turbulent spectral 
slope, and any correlations between components of v. 
As explained in ^A.2l j is distributed as a Maxwellian, 
or, crudely, 



logio fj = 



-0.16 
-0.49- 



(11) 



The above results pertain to velocity fields with 
strictly uncorrelated components, which is not en- 
tirely realistic. An uncorrelated flow can be decom- 
posed into an irrotational, potential, compressive por- 
tion and a vortical, incompressive portion. In Fourier 
space, potential flows are aligned with the wavevector 
and vortical flows are perpendicular to it. Since there 
are twice as many degrees of freedom in the perpen- 
dicular direction, an uncorrelated flow contains two- 
thirds of its power in vortical motion and one-third 
in potential flow. Since core angular momentum de- 
rives entirely from the vortical part of the flow, it 
is Y^3/2 times higher in incompressible turbulence 
as compared to uncorrelated turbulence of the same 
a. Equation Q implies turbulent Mach numbers of 
about unity on the core scale, suggesting behavior in- 
termediate between the incompressible and pressure- 
less limits. It is reasonable, therefore, to multiply fj 
by 1.11 (the harmonic mean of ^3/2 and unity) to 
account for correlations introduced by pressure. This 
correction is however rather minor. 

2.2. Collapse of Cores 

Cores collapse with the accretion rate ch aracteristi c 
of any self-gravitating body, M = c^g/G (|Shdll977f ). 
We assume only a fraction e of the core successfully 
accretes, so the final stellar mass is M = eMc, and 



£ 



G 

2.8 X 10"^ 



f 



3/4 



p3/2 
■10 



yr 



(12) 



There exists an initial phase of more rapid accre- 
tion during which the central region of the core col- 
lapses homologously; however, equation (fT^ is valid 
for most of the mass. 

The theorv p redicting e ~ 25% — 80% 
( Matzner fc McKeei.2000. ^ invokes the removal of ma- 
terial from the rotational axis by protostellar jets; 
this somewhat increases the mean angular momen- 
tum of accreting material, but we make no correction 
for this effect. 



It is important to note that equation (|12() relies 
on the assumption that cores are the mass reser- 
voir for stars and collapse in their entirety. This 
accounts neither for the possibility that accretion 
migh t be aborted before the collapse has finished 
(e.g.. lEeipurth fc ClarkdEoOll ). nor for the possibil- 
ity that accretion continues after the core has col- 
lapsed (Bonneh et al. 200 J). These effects would al- 
ter the amount, duration, and angular momentum of 
accreted material; however, we assume they can be 
neglected when estimating the characteristic scales 
of star formation. 

The characteristic disk radius formed during infall 
is f/{GM), or 



/I 



Abef; 



1/2 



Rn 



,3/2 

Jp 



^cl 



(13) 

This is comparable to the critical radius beyond 
which the disk's self-gravity can cause local fragmen- 
tation, at least in principle; see equation H31|) below. 
The maximum disk period is of order 9000 years, but 
is quite sensitive to the parameters fj, e, and fp. 

The inverse relation between R^ and e results from 
the fact that a decrease in e does not reduce the spe- 
cific angular m omentum j of accreted ma terial. This 
is true for the Matzner &: McKe3 (|20o3) model, in 
which an outflow removes material only along the ro- 
tational axis. In scenarios where e < 1 because accre- 
tion is aborted, the disk inherits only low-j material 
from the core center and R^ (and the tendency to 
fragment) are accordingly suppressed. 

2.3. Comparison to Prior Models 

We have not attempted to model the internal struc- 
tures of cores under the influence of both thermal 
and nonthermal support, as previous authors have 
done (iMvers fc Fulleiill992l : lMcKee fc Hollimanlll99fll : 
Currv fc McKeelll999^ . ' Our emphasis has instead 
been on the overall mass, radius, and angular mo- 
mentum scales of the smallest collapsible objects that 
can exist in a given turbulent environment. For this 
we rely on a Bonnor-Ebert model, rescaled to account 
for nonthermal support, and neglect the flattening of 
the density proflle this nonthermal pressure should 
provide. Although far from perfect, this approxima- 
tion is best for the most thermally-supported cores 
on which we concentrate. 

In our studv of core rotat ion we have followed 
Burkert Bodenheimeil 12000) in positing a homoge- 
neous turbulent velocity field obeying the line width- 
size relation ((2| that pervades cores and gives them 
rotation. Studies of core internal structure have of- 
ten assumed instead that the nonthermal pressure is 
a function of the local density, e.g., in constructing 



6 



multi-pressure polvtro pe (iMcKee &: Holliman 
or nested polytrope ( Currv &: McKed Il99^ ) core 
models, or that it is a function of the local d istance 
from the core 's center (|Mvers fc Fulled Il992ll . Poly- 
tropic models are motivated, in part, by the poly- 
tropic behavior of Alfven wave pressure in the WKB 
approximation ( McKee fc Zweibelll995l ). This is an 
idealization because of the predominance of long- 
wavelength p erturbations and because nurnerica l sim- 
ulations Ce.g.. lStone et a1.lll99Sl : lMac T.owl ll999h indi- 
cate that Alfven waves damp rapidly in molecular 
cloud conditions. Our assumption of homogeneous 
turbulence is equally idealized: it does not explain 
the origin of the turbulent field, and it ignores the 
existence of a line width-density relation. We also 
have not accounted for the possibility of a break in 
the line width-size relation on scales between the core 
and its parent clump or cloud. Nevertheless, our re- 
sults are in reasonable agreement with observations 

Our model is similar in some ways to the 
TN T (thermal plus nonth ermal) model introduced 
bv iMv ers fc Fuller fl992!) and used recently by 
McKee Rr. Tan (2Q03.). However the TNT model 
posits a nonthermal line width that depends on the lo- 
cal radius from the core center: it is akin to the poly- 
tropic models in that it associates turbulence with 
the local conditions. The TNT model also idealizes 
core density structures as singular power laws, rather 
than accounting for the flat central region that exists 
prior to collapse: this is significant because it implies 
that thermal pressure will always dominate within 
some radius. In contrast, our models argue for com- 
parable thermal and nothermal contributions even in 
maximally thermal cores. 

Our predictions for the mass and radius of precol- 
lapse cores, equations ©, ar e similar to tho s e of t he 
smallest (thermal) cores in iMcKee faril (j2Qo3)'s 
theory. However these authors do not account for 
the possibility of a minimum level of turbulent sup- 
port. As a result, their thermal cores are less mas- 
sive and more slowlv accreting th an we would p re- 
dict. Although Mvers fc Fuller and McKee fc Tanl do 
not address the core angular momentum distribution, 
their thermal cores would be significantly more quies- 
cent and more slowly rotating than the cores modeled 
here. 

2.4. Comparison to Observations 

Masses: The mass and radius scales derived in ^ 
(eq. ini) are consistent with the typical properties 
of prestellar cores obs erved in near bv star- f orming 
regions, like Orion B (|,Tohnstone et al ] l2nnih . The 



power law to masses above the minimum thermally- 
supported mass identified in this section. 

Nonthermal linewidths: Ijiiina. Mvers. &: Adams! 



()1999I ) classified presetellar cores according to their 
association with cores and IRAS sources. In their 
database, those cores that are truly starless in this 
comparison have median ratios of nonthermal to ther- 
mal linewidths about 35% smaller than suggested by 
equation (jHI - see their figure 4. The observed distri- 
bution overlaps our prediction. Some of the discrep- 
ancy may be attributable to rotational motions that 
do not contribute to the (beam-resolved) linewidth, 
o r to static magnetic fields . 
Jiiina. Mvers. &: Adamsl observe a correlation be- 



tween the thermal and nonthermal linewidths: o"nt oc 
a^^^ (their figure 7 and table B7). At the cold 
and quiesent end of this corre lation are L1 4 98 an d 
L1517B, studied in detail by iTafalla et all ((200i). 
This trend is not explained in our theory, in which 
equation Q gives unt — Cth- 

Recall however that equation Q only represents 
the characteristic scales of cores. It relies on several 
assumptions: 1. that the core is at the brink of col- 
lapse; 2. that it is free of magnetic support; and 3. 
that it has the minimum level of turbulence consis- 
tent with assumptions 1 and 2 in the context of a 
homogeneous model for the turbulent velocity field. 
An individual core need not adhere to these assump- 
tions, although we have argued that they suffice to set 
characteristic scales. In the case of the very quiescent 
cores^ L1498 and L1517B, we suspect that magnetic 
fields rather than turbulence are supplying a signifi- 
cant amount of nonthermal supp ort. 
Nonthermal support: ..Johnstone et al.l ( 200lh identify 
an effective temperature, using Bonnor-Ebert, mod- 
els which is typically 20-40 K. Many of these clumps' 
actual temperatures measured at less than 20 K, with 
10 K being typical. A significant contribution of non- 
thermal pressure - as in our equation Q - offers an 
explanation for this discrepancy. 

It is interesting to note that this level of nonther- 
mal support holds even in the very quiescen t cores 
L1498 and L1571B, for which iTafalla. et ^^\\ ()2n04l ) 
find Teff = 28 and 20 K, respectively, in their Bonnor- 
Ebert fits. Compared to the measured thermal tem- 
peratures Tc = 10 and 9.5 K, we note that the ratio is 
consistent with equation Q although the linewidths 
are not. Although part of the nonthermal support of 
these cores is rotational, most is likely magnetic. 
Rotation rates: As for core rotation rates, equation 
(fTUl) agrees with th e value j ~ 10^^ cm^ s~^ deduced 
bv lGoodman et al. ( 19931 ) f or cores of similar mass. 
The s econd line differs from lBurkert fc Bodenheimed 
(|200nl )'s formula j = 7 x 10'^'^{Rj0.1pcf/^cm'^ s'^. 



core mass spectrum in such surveys typically shows a 

^These are not characteristic of the 172 lJiiina. Mvers. fc Adamj lUgogi) cores: for instance, L1498 is their smallest core 
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However, Burkert &: Bodenheimeit s results were nu- 
merical; our own unpublished numerical experiments, 
which were based on their method, exhibit a bias to 
underestimate j by about the same factor. We believe 
this arises due to the periodic boundary conditions 
employed in the numerical grid; due to limited dy- 
namical range, these impose correlations on the box 
scale that affect the spectrum as sampled on the core 
scale and smaller. 

The rotation rates observed bv lTafalla e HIl ^2004 ) 
for the very quiescent cores L1498 and L1517B are 
only a factor of a few smaller than the typical val- 
ues predicted by equation (fTU)) . and are within the 
distribution descri bed by equation (ITTl). 
Accretion rates: iMotovama &: Yoshidal ( 20031 ) re- 
view a range of observational evidence that mass 
accretion can significantly exceed the characteristic 
value Cg/G in an early (class 0) accretion phase. 
They suggest the initial phase of homologous col- 
lapse, or collapse induced by an external impulse 
(see also iHen nebelle et al. 200^), as possible expla- 
nations. These possibilities are both very reason- 
able, but note that we predict an accretion rate of 
roughly 2.6cf/G, on account of the nonthermal sup- 
port (eq. ^2]). Therefore, the need to invoke alter- 
native explanations for the observed protostellar ac- 
cretion rates is partially alleviated. 
Magnetic fields: Finally, we comment on the 
neglect of static magnetic fields in our calcula- 
tions. Employing the Chandrasekhar-Fermi method, 



mg tne 

Crutcher et al.1 (|2004l ') find that static magnetic fields 
in prestellar cores could either be subdominant or 
comparable in magnitude with thermal pressure, de- 
pending on the application of a geometrical correc- 
tion factor to their observations. Theoretically, one 
expects that ambipolar diffusion will sap the influence 
of magnetic fields prior to core collapse. In addition, 
infere nces of short core lifetimes (e.g., IVisser et al 



l2002h and obs e rvatio ns of high core column den- 
sities ( Nakanol ^^3) sue taken to indicate that 
static magnetic fields are not significant in core sup- 
port. On theoretical grounds, Zweibel fc .TosafatssonI 
find that ambipolar diffusion is accelerated 
in turbulent gas. We conclude that our approxi- 
mation of zero mean field is adequate given current 
uncertainties, while recognizing that this should be 
re-examined in future work. As discussed above, 
magnetic fields may be especially important for un- 
derstanding the thermal- nonthermal l inewidth cor- 
relation seen bv lJiiina. Mvers. Adams (1999 ) and 
the nontherm al sup port of quiescent cores seen by 
Tafalla et al.1 (|2004l ). 



The first opportunity for fragmentation is in the ini- 
tial phase of collapse, which follows from the homol- 
ogous contraction of the core's central region. In this 
phase, the mass accretion rate can exceed equation 
(|12j) . If the core's central region rotates sufficiently 
rapidly, it will form a disk that can fragment. 

This phase of fragmentation has been studied in 
detail recently by Matsumoto & Hanawa (20qJ), who 
consider the ini tial collapse of slowlv ro tating Bonnor- 
Ebert spheres. iMatsumoto fc Hanawal find that frag- 
mentation occurs by one mode or another so long as 
^cts,ctT ^ 0.045, where ilc is the rotation rate and 
%,ctr is the free-fall time evaluated at the middle of 
the initial state. Using equations ®, (fTU)) . and the 
moment of inertia Ir 



0.2827Mci?2, we find 



^^c%,ctr - 0.25 



Jp 



(14) 



By this criterion, then, all but the most slowly rotat- 
ing cores (/j < 0.18/p''^; about 3% of the population 
according to eq. |A20j ) should undergo fragmenta- 
tion in the initial collapse phase. 

Ho wever, this conclusion i s not r obust. Note, first, 
that Matsumoto &: Hanawal (20o3) rely on the ap- 
proximation of barotropic core collapse (isothermal, 
for nf{ < 10^^ cm~^) rather than a det ailed thermo- 
dynamic calculation. IBoss et aP (|2000l ) caution that 
protostellar collapse simulations employing the Ed- 
dington approximation differ significantly with regard 
to fragmentat ion from thos e that use a barotropic 
law. Recentlv. iLarso P (|2004l ) has argued that thermal 
evolution plays a critical role in determining stellar 
masses. Likewise, nu merical arid ana l ytical studies 
of disk fragmentation (Gammie I 2OOII : iPickett et al 



[2003, : ,Rafikov.2004 : ■Levin.2003. ) stress the importance 
of cooling. We will return to this issue in detail in ^ 
and add some comments in ^6.21 



3. PROMPT FRAGMENTATION 



Second, IMatsumoto Sz Hanawa 's calculations stop 
after at most O.O7M0 of several Mq have accreted 
into fragments. As they point out, the future evo- 
lution of these fragments is difficult to predict. We 
note that the remaining mass accretes through a disk 
around the newly formed objects. The specific angu- 
lar momentum of this material greatl y exceeds that of 
the ob served fragments. As shown bv Ba te"&: Bonnell 
(|l997l ). a binary tightens due to gravitational inter- 
action with the circumbinary disc. We suspect that 
this drives prompt protostellar fragments to merge. 
Future simulations will be required to test this hy- 
pothesis, especially for the complicated case of multi- 
ple prompt fragmentations. If prompt fragmentation 
does succeed, it should produce tight binaries with j 
significantly below the typical value cited in equation 
()10() - an therefore, significantly below the median of 
the binary distribution ((Duaucnnov &: Mavor.,199'll ). 
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For both of these reasons we take equation ((11]) to 
be only suggestive of fragmentation in the initial, im- 
pulsive collapse phase. We probe fragmentation in 
later stages of accretion in the upcoming sections. 

4. STEADY PROTOSTELLAR DISKS: LOCAL 
GRAVITATIONAL INSTABILITY AND LOCAL 
FRAGMENTATION 

We now explore the stability of the steady proto- 
stellar accretion disks that form during the collapse 
of cores like those described in ^ during the phase 
of main accretion during which the accretion rate 
is given by (|12() . There are two modes by which 
self-gravitation can trigger disk fragmentation, one 
of which is su bject to a local instabilit y criterion 
(lloomrel and the other a global one (!Sh u et al.l 

1199(1 ). We focus initially on the local instability and 
return to the global one in ^ The consequences of 
Toomre's instability have be en greatly elucidated by 
the numerical simulations of iGammid (200l|), as de- 
scribed below. 

iLin fc Pringi3 (1199(1 ) have carried out an analysis 
of protostellar disk evolution, based on initial con- 
ditions qualitatively similar to those sketched in ^ 
They implemented a gravitational viscosity that rose 
as Toomre's instability parameter. 



vrGS' 



(15) 



declined through its critical value of unity. Here O 
and S are the disk's orbital frequency and column 
density, respectively; and 



ldJ_ZPdz 



5S 



the other hand, strongl y unstable disk s enter a state 
of "gravitoturbulence" (lGammiell2nnih in which the 
coherence length is only a few effective scale heights 



In this state angular momentum tra.nspo rt is a lo- 
cal process (see also iLodato fc Ricd l2004^ . At the 
brink of fragmentation. IShakura fc Sunvaevl's viscos- 
ity pa rameter a takes maximal value, which ICammid 
(WM) shows to be^ 



a 



9(7; 



2d 



i)nt 



cool 



(16) 



if ^cooi is the cooling timescale in the critical state. 

We recall that a is related to the mass accretion 
rate through the disk: 



(17) 



The subscript "vise" is added to distinguish the ac- 
cretion rate due to local viscous transport from, for 
i nstance, the rate of mass supply (eq. ^2] ) • 
GoodmanI (120031 ) points out that angular momen- 



(where z is altitude) is the signal speed in the disk, 
which is slightly above the midplane isothermal sound 
speed Cs = yJPj p. With this prescription for an- 
gular momentum transport, they found Q saturat- 
ing typically in the range 0.2-0.5. This study lacked, 
however, two physical important physical effects: 1. 
the nonlocal nature of gravitational angular momen- 
tum transport when Q > 1; and 2. the dissolution 
into fragments of disks with Q <^\. Both shortcom- 
ings are addre s sed b y th e nonlinear sim ulations of 
Laughlin etlll ((l99i) and'Ga mmi3 (l200lh : 

1. When gravitational instability is weak and Q > 
1, spiral modes saturate at low amplit udes via mode- 
mode coupling ( Laughlin et al.iri99 <^: in this state, 
redistribution of angular momentum is nonlocal and 
is described onlv approximately bv a local prescrip- 
tion, i.e., by IShakura fc Sunvaevl 's a parameter. On 

''We define a using the isothermal sound speed Cs, rather than 
differs by a factor 72d relative to its definition in iGammiei (.200 . 



tum loss in magnetized winds can lead mass to accrete 
more rapidly than equation ((T^ predicts. We how- 
ever assume that strong winds do not develop over 
most of the disk radius, and we adopt equation ((TB|) 
maximal value. 

2. By identifying the precise boundary between 
gravit oturbulent and fragmenting disks, iGammij 
(|200ll ) provides the means to discriminate stable mass 
accretion from the formation of new bound objects. 
A number of papers argue that the precise outcome 
of fragmentation depends s ensitively on the th ermal 
properties of the disk fe.g.. IPickett et aP 200,"^ ): nev- 
ertheless the threshold for instability is given ade- 
quately by equation ((20|) . 

These ins ights reveal two points about 
iLin &: Pringld 's simulations. First, they estimate in a 
reasonable fashion the dynamics of protostellar disks 
in their main ac cretion ph ase (since their estimate of 



a approximates Gammid 's result when (5 = 1); and 
second, that these disks are close enough to Q = 1 to 
merit an careful investigation of whether fragmenta- 
tion does indeed occur. 

The results of Gammie's two-dimensioi ial simula- 
tions (corroborated in three dimensions bv lRice et al 



'2003^ provide a practical, local description for the 
disks that have entered a state of self-gravitating tur- 
bulence. This state separates disks that are relatively 
smooth from those that have fragmented into swarms 
the disk signal speed Cd, in equation 1171 . As a result our a 
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of clumps; it also is a state in which the mean mid- 
plane conditions (density, temperature, etc.) are rea- 
sonably constrained. We use it as a litmus test for 
fragmentation in a variety of circumstances. 



4.1. Local criteria for fragmentation by Toomre's 
instability 

Toomrel (Il964 ^ demonstrated the onset of axisym- 



metric instability in disks with Q < 1 ; although other 
modes grow for somewhat higher O. iGammid ()20nil ) 
verified that Q ~ 1 at the boundary of fragmentation. 

To assess whether the disk can process accretion 
without fragmenting, we adopt 



l2dCs 



where 



l2d 



91nS 



(18) 



in equat ions (11511 and (11611. In the limits of weak self- 
gravitv jGoldreich et al.lll986^ or strong self-gravity 



( Gammiell2nniir repsectively. 



72d 




(Q » 1) 
{Q « 1) 



(19) 



w here 7 is t he ad iabatic index of disk gas. 

Gammiel ()2flnih adopted the specific value ■j2d = 2, 
corresponding to 7 = 3 or 2 for strong or weak self- 
gravity, respectively. Either value is stiffer than a 
realistic disk. In the temperature range of interest, 
a protostellar disk obeys 7 = 5/3 because rotational 
degrees of freedom are not excited;^ then 72^ = 3/2 
or 9/5 in the strong and weak self-gravity limits, re- 
spectively. At the onset of fragmentation Q ~ 1, and 
we assume that 



l2d - 1-65, 



(Q = i) 



an intermediate value. We also assume that the crit- 
ical cooling time remains roughly one half orbital pe- 
riod, i.e., r^fcooi = 3, for the softer equation of state; 
this remains to be checked by future numerical ex- 
periments. With these choices, equation (fT6|) gives 



a 



0.23. 



(Q = l) 



(20) 



Weak and strong self-gravity give a = 0.30 and 0.19, 
respectively. We adopt a = 0.23 clS cl fiducial maxi- 
mal value, but consider this uncertain by about 50%. 

The condition Q = 1 may immediately be ex- 
pressed in terms of c^, the midplane density p, and 
the midplane pressure p: 



1/2^' 

l2d ^ 



l2d 



27rG 



and p 



TT 



272^1/2 



(21) 



where the middle expression uses the approximate re- 
lation p = Sr2/(2cs), and the last expression uses 
p = p(?g. These relations are only approximate, 
as they neglect the disk's self-gravity in its vertical 
structure, but they suffice to delineate the onset of 
instability. 

An important additional criterion comes from com- 
bining eqs. ((T^ . (|T3|l. and ((T7|l to give 



Q 1/2 



eQ 



(22) 



-cff 



Conversely, if the disk processes gas in steady state 
(M^isc = Mace) then 



Q 



(steady — state) 



(23) 



Using 72rf = 1.65 and a = 0.23 to describe the onset 
of fragmentation, we find that the disk can process 
steady accretion with Q > 1 only when 



(24) 



i.e., using eq. (jHI to describe a typical unstable core. 



T > LOSe^/^Teff ^^(-^, 



fl'^ Uo%j 



\ 2/3 



(25) 



(where T refers to the midplane temperature); oth- 
erwise fragmentation will occur where the gas falls. 
For typical core temperatures of order 10 K, disks 
colder than about 15.6 K are unable to process the 
high accretion rate and will fragment into clumps. 
This critical temperature is rather insensitive to the 
uncertainties described after equation (|20(1 : it varies 



as a 



4.1.1. Self-luminous disks 

If the heat generated locally by viscosity is respon- 
sible for the midplane temperature, then criterion 
(Pljl may be evaluated di rectl y. Our treatm ent fol- 
lows E^Priniii 11990') and iLevinI (j2003l 'l. among 
others. In thermal steady state, the flux of viscous 
energy radiated by each face of the disk is related to 
the accretion rate through 



*Shock heating and radiative transfer alter the effective 7 from 
However it is the adiabatic 7 that enters in the theory of lGarnrnia 



3J72M 
Svr 



(26) 



its adiabatic value, especially for disks close to fragmentation. 
.2001.1. 
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The flux can also be derived from radiation transfer 
across an optical depth from the disk midplane 

to its surface; therefore 



aT^ X 



4tpi (t < 1) 
8/(3rR), (t»1) 



(27) 



Q ^ 1 (eq. ^ni)- We find that fragmentation ensues 
when 



Cs < Cs,T (r <. I) ^3Q^ 
Cs > Cs,r (t > 1) 



where tr,pi = kr_piS/2 is the optical depth corre- 
sponding to Rosseland or Planck opacity, or Kpi, 
respectively.^ The factor 16/3 arises in the optically 
thick case when t he dissipation rate per unit mass is 
a constant (e.g., Chick &: CassenI and is only 

approximate. For the disk to achieve an equilibrium 
temperature in which local heating balances radiative 
cooling, we must have Fy = F^. 

The midplane temperature and sound speed are re- 
lated via 

4 = ^ (28) 

where fi ~ 2.3mp is the mean molecular weight. We 
neglect the radiation pressure in this expression. This 
is entirely justified for present-day star formation, 
since when Q = 1 the ratio of gas to radiation pres- 
sure is 




4.3 X 10 
JkR 



-10 



fi /0.23 1.65\V3 _^ 



2.3m.o 



(31) 



corresponding to a period of 460 years for the fiducia l 
parameters (400 yea rs for [Alexander Sz Fergusoiil 's 
opacity; 350 years for Henning Sz Stoenienkol 's opac- 
ity), and 



0.174 



j.l/30 j.1/10 
JkR JkPI 



A* 



a 
023 

8/15 



1/15 



1.65 

l2d 



1/30 



km s 



(32) 



3k B^'^ 
27rGa/ir3 



/lOOOyr 
V period 



500 K 



using eq. (ETJ. 
For the opacity we adopt 



3.0/«R 
8.4/«pi 



X 10"*r^ cm^ 



This describes the ice branch of Semenov et al 



129) 



()2003l Vs ice grain opacities in their "composite ag- 
gregate" model; ea. (I29II fits their results around 



16 K. Note that [Alexander fc F erguson (1994) com- 
puted somewhat l ower opacities (/kr ^ 0;67; see also 
Bell fc Linlll99i 'l. iHenning fc Stognienkd also 



To summarize, disks with opacity law (|29() that are 
warmed solely by their own viscosity will be subject 
to local fragmentation in any region whose orbital 
frequency is less than the value given in eq. (|3()|) . 

It is worthwhile to emphasize two points about the 
above equations. First, Qcr is an upper bound for the 
frequency of a locally unstable disk, i.e., all locally 
unstable disks have periods of roughly 460 years or 
longer. This bound depends only on the opacity law 
of the disk midplane and is independent of the mass 
supply rate and disk surface density (so long as the 
disk is optically thick). This result is specific to a 
cancellation that occurs when the opacity varies as 

Second, the lower bound on the period at which 
fragmentation can occur is pushed to even higher val- 
ues if the disk is optically thin or if it is affected by 
external irradiation. The division between optically 
thin and thick fragmentation depends only on the 
midplane sound speed c^; the critical value Cg^r cor- 
responds to a midplane temperature of roughly 8.4 
K. But recall that the fragmentation temperature is 
determined by the mass supply rate to be roughly 
T ~ 16 K for cores with = 10 K (eq. j25]V 
This implies that disks are typically opticall thick 

^Recently. iJohnson fc Gammid ll2003fl have published numerical experiments on razor-thin discs and have shown that Eq. 1271 
is generally inaccurate when the disc develops strong gravitationally-driven turbulence. However, in their model the breakdown of 
Eq. 1271 is serious when the opacity is a sensitive function of temperature, which occurs only when the dust begins to evaporate 
at above 1000 K. Temperatures this high are clearly irrelevant for protostellar discs beyond a fraction of an AU, and thus for the 
pro blem at hand Eq. 12711 is adequate. We also note that when the opacity is temperature sensitive, then the one-zone model of 
iJohnson fc Gammia does not treat faithfully the vertical radiative transport in a real disc. 



derive a lower value (/^r ^ 0.43) after allowing for 
dust agglomeration in the parent molecular core. We 
include the two factors to account for variations 
in metallicity, dust processing, and uncertainties in 
opacity modeling. 

Equations (HH), (|17|, and ^ provide enough 
information to solve Fy = F^ for T(r), given an accre- 
tion rate M. As discussed above, the local accretion 
rate (eq. ^7j) must match the supply (eq. ^2). To 
determine the onset of fragmentation, we also impose 
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at the fragmentation bound Indeed, using equations 
(Uni) and (jSH), all cores with 



Te > 5.4 



A/5 .1/15 
/ftPl JkR 



K 



(33) 



will produce disks that are significantly more prone to 
fragmentation than those from colder cores because 
of the transition to optically thick disks. (In equation 

we suppress the variation Tc oc a^/^^2d^^ f^~^^^^ 
around fiducial values of these parameters.) 

Equation (|30j) argues strongly against the forma- 
tion of giant planets by direct disk f ragmenta t ion, i n 
agreement with the recent work by iRafikovl ((2QQ^ . 
This argument gains strength when the effect of stel- 
lar irradiation is considered. 

Fragment masses. We pause here to note the char- 
acteristic masses of fragments formed by local gravi- 
ta tional instabi lity. This topic is investigated in detail 



tationai mstaDiiity. i nis topic is investigated , m acta: 
bv lLevinI (j2003) and lGoodman k Tanl(l2nn4^ . If frag 
ments form, their initial mass is of order S(27r/f)^, 
i.e., roughly Jupiter's mass for both the optically 
thick and optically thin cases if Tc = 10 K (but scaling 



as Tc^'^ and Tc ''^ respectively). Continued growth 
is likely up to the gap-opening mass, which we es- 
timate at 36 Jupiter masses for fiducial parameters, 
and perhaps beyond. Local fragmentation is t here- 
fore most relevant to the origin of planetary ( Bosd 
1993) or brown dwarf fiBate et al.i r2nn2) companions. 



-7/2 



4.2. Irradiation 

The preceding section identified a critical midplane 
temperature for local fragmentation (typically ~ 15.6 
K in disks around low-mass protostars, eq. [2S1), and 
a critical disk period (typ. ~ 460 yr) above which vis- 
cous heating permits cooler midplane temperatures. 
However as we mentioned before, the argument in the 
preceding section neglects external irradiation of the 
disc. Here we self-consistently take irradiation into 
account and show that it completely prevents local 
fragmentation for a wide range of core parameters. 

To illustrate its importance, suppose irradiation 
normal to the disk surface is a fraction fp of the 
spherical flux. The equilibrium temperature of an op- 
tically thick disk is then given by crT^ = /irL/(47ri?^), 
if L is the accretion luminosity; or 



94 f^/^ f^y^' /460yry/^ 
JF \ AKnn I ^period/ 



4500 
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thin, its temperature is generically higher than given 
above: for a local radiation energy density is aT^^ 
and color temperature is Tc, the equilibrium temper- 
ature T of an optically thin region is 



T^KPi(T) = Tgt Kpi(Tc 



(35) 



where Kpi refers to the Planck average opacity. Since 
Kpi(T) is monotonic in the temperature range of in- 
terest, Tegy < T < Tc. Furthermore, since the fiux 
at the disk surface cannot exceed aT^^yC, an optically 
thin disk must be warmer than a thick disk in the 
same radiation environment.^ 

The fiux suppression factor fp can be quite 
small in the case that the disk is illumi- 
nated at a grazing incidence from the stel- 
lar photosphere (e-g.. iKenvon fc HartmannI Il987l : 
Chiang fc GoldreichI 1997^ . For an accreting pro- 
tostar, however, reradiation from the stellar en- 
velope tends to illuminate and heat the disk 
Keene k Ma,sso'^ 1199(1: iNattal IT99^: IChick et a] 



1 99Bc IChick Cassenlh 9971: iD'Alessio et aklll 99' 

To gauge this, consider the column of the infall en- 
velope. The spherically averaged column density out- 
ward from some radius r is Esph(r) = M/[27rr7;esc(r)] 
Evaluated at the centrifugal radius in eq. ((T^ . 



5^sph(^d 



0.43 



.3/4 
Jp . 



(36) 



2 X 10- 



'M0yr-iy 
(34) 

where the stellar eff'ective temperature is T^. (Note 
this is independent of the stellar mass M, if ref- 
erenced to orbital period.) If the disk is optically 
^The lower limit on T in the thin case is y/2 below that in eq. 



In spherical infall Esph(r) diverges inward as r~^/^. 
For rotating infall, however, only low-j material per- 
sists within i?^; this makes Ssph(i?(i) a characteristic 
value for the column inward from R^, as well. 

Intriguingly, the characteristic column T,sph{Rd) is 
roughly Sci/2, i.e., half that of the parent clump 
or cloud, at the end of accretion. (Beforehand 
it is much greater: the column to the disk edge 
varies as t~^ during steady accretion.) Deviations 
from spherical symmetry due to angular momen- 
tum are given by the rotating-infall solution of 
iTerebev. Shu, fc Casseiil (Il984l 'l. Additionally, the 
protostellar jet will scour the axis to produce a cavity. 

Since star formation appears t o be limited to re- 
gions of hi gh visual ext inction ( Onishi et al. 19981 : 
iHara et al.lll99 9: Johns tone et al.l 2004, as suggested 
theoretically bv (McKeei 1989), and since T,sph{Rd) — 
Sci/2, starlight is absorbed and reprocessed after it 
traverses a fraction of the infall column. This repro- 
cessing is strongly aff'ected by the presence or absence 
of an outfiow cavity: therefore, we consider both cases 
below. 

13411 ■ but this factor is not significant. 
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The temperature of an irradiated disk depends on 
the optical depth of the disk to the incoming radia- 
tion and to its own thermal emission. For a disk at 
the fragmentation threshold with Q = 1, 



3q /if 



23/2 



210 



A/2 , 
Jp / 



ff V60% 



4/3 



(37) 



The first hne derives from equations (fT5|) . (|T8|) . (|2i|) . 
and (|36j) ; the second line applies the results of ^ (and 
varies as (72d/a)^^^ around our fiducial values). The 
outer disk column is thus of order 90 Sd, or ~ 3 g 
^ for typical clouds. According to the opacity 



cm 



law ()29() . the disk will be vertically optically thick to 
Planck radiation of color temperature exceeding 



19.6 



f 

■'3 



3/2 



j.1/2 „5/8 
JkP\ h 



( 165M0pc-2~ 



1/2 



K. 



The disk's opacity to radiation from dust of a cer- 
tain temperature is higher than to a Planck field of 
that color temperature: /^pi — 2.5 for dust radiation. 
This lowers the critical illuminating temperature to 
about 12 K, and oblique illumination lowers it further 
still. The disk can therefore be assumed optically 
thick to irradiation, although it is only marginally 
thick with respect to its own radiation at a fragmen- 
tation temperature of ~ 15.6 K. 

4.2.1. No outflow cavity 

The reprocessing radius typically lies well inside Rd 
but outside the dust evaporation radius, because the 
latter is optically thin in the end phases of accretion. 
Disk irradiation is therefore due primarily to the sec- 
ond reprocessing of starlight at radii around Rd- We 
shall identify the dust temperature at the reprocess- 
ing surface, which affects the amount of re-absorption 
and re-emission above the disk. 

Close to the pro tostar, the angle-averag e d col - 
umn density in the iTerebev. Shu. &: CassenI (|l984l ) 
infall profile is S(r) = 1.17 {r / R^Y I '^T.^^Y.iRd) + 
0{r/Rdf'^. We identify the reprocessing radius 
ri with the surface of optical depth unity, then 
S(ri)K^ = 1, where = K-p\{Ti,) is the opacity to 
starlight of color temperature T^. In addition, equa- 
tion (|35|) determines the dust temperature Ti at this 



location: 
these, 



Kpi(ri)Tf = 



K,^L/{4:Trrfac). Combining 



' 26.7 

''Ti also depends weakly on the gas density at ri, 
ignored in eq. EH 



(38) 



The left hand side is a somewhat complicated func- 
tion of Ti because the dust composition, a function 
of Ti, determines k^, and Kpi(ri); metallicity and the 
stellar temperature T*. also enter. ^ The right hand 
side depends on model parameters. A self-consistent 
solution of (|38|) is not always available; this signals 
a change in the dust composition and fixes Ti at a 
composition boundary. Nevertheless, this equation 
demonstrates that the reprocessing temperature Ti 
is a pure function of LT,'^^^ for fixed T*. and metallic- 
ity, so long as the infall is optically thick to starlight. 
In a range of Ti for which the dust composition is 
constant, Ti oc L^/^Sgpj| because Kpi(Ti) oc T^. 

Given a optically thin emission at the reprocessing 
temperature Ti, the effective opacity for the second 
reprocessing event is given by 



^ K^iB^{Ti) 



(39) 



where Kyi is the opacity law at the reprocessing ra- 
dius and is that for cool dust near Rd- Figure ^ 
shows Ti and R for a given model; notably, k ~ 15 
cm^ g~^ for a wide range of parameters, dropping 
only for especially large disk radii (and especially low 
values of T.^p^{Rd)). 

Integrating the flux through the disk at Rd from 
the sky intensity of reradiated light (assuming opti- 
cally thin second reprocessing), we find that the flux 
of this radiation at the disk radius Rd is given by 
equation (jHH) if 



fp = 0.28KEsph 



(40) 



giving fp ~ 1/16 for typical parameters. Using this 
result to calculate the irradiation temperature of an 
optically thick disk in eq. ()34() . we find that it falls be- 
low the critical temperature for local fragmentation 
at the outer disk edge only if 



period > 



7400/pi 



1.69 



3/4 



^cl 



V4500K 



10cm2 g-i 165M0PC-2 
10K\9/4 



yr 



3/4 



(41) 



there; this is comparable to the characteristic maxi- 
mum period in H2.2L 

Note we have taken e = 100% here, since the ab- 
sence of the outflow cavity implies that all of the ma- 
terial is accreted onto the star. This raises the critical 
disk temperature to about 22 K. 
which afltects the evaporation temperatures of the dust components; this is 
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Fig. 1. Lower panel: Temperature Ti of the reprocessing s urface in envelopes lac king outflow cavities. Ti is computed as 
a function of outer disk period period using the opacity model of lSemenov et alT (|200.?t) . Assumed parameters are M = O.SMq, 



Af = 3 X 10~^Mq yr~^, T* — 4500 K, and solar metallicity; however, equation 1)38^ permits this diagram to be rescaled for other 
parameters. (The density dependence of evaporation temperatures has been suppressed for simplicity; lines here are drawn for 
10~^^ g cm~^.) Upper panel: effective opacity k of cool dust to optically thin thermal radiation from dust of temperature T\. 



4.2.2. Outflow cavity 

A more realistic calculation must account for the 
region along the axis that is cleared by the jet. An 
outflow cavity warms the disk by providing more di- 
rect illumination than is available otherwise. To cal- 
culate this we must specify the shape of the outflow 
cavity. 

Note, first, that the wind ram pressure is expected 
generically to vary as 1/fr^ sin^ 9) with distanc e and 
angle from the axis (jMatzner McKeelll 99fll V In- 
flow ram pressure scales as r~^/^ from the centrifu- 
gal radius to the edge of the inflow. We see on 
comparison that inflow dominates close to the star 
and near the equator, and the wind dominates far 
from the star and near t he axis. By this reasoning, 
Matzner fc McKed (1200(1 ) divided the initial core into 
accreting and ejected angles depending on the veloc- 
ity imparted by the wind impulse to gas at the edge 
of the core. Our parameter e is simply the accreted 
mass as a fraction of the initial mass. 

By the same logic, gas that is not cast away by the 
wind is destined to fall inward, and its motion is less 
and less affected by the wind ram pressure as it does. 
Therefore we are justified in approximating the shape 
of the outflow cavity as an unperturbed streamline i n 
the infall solution of lTerebev. Shu, fc CassenI 
The streamline in question is roughly the one that di- 
vides infall and outflow in Ma.tzner & McKee ( 2000)'s 
theory. For the case of a spherical initial core, then, 
the initial angle of this streamline is given by 



(We shall assume a spherical core for the remainder 
of this section.) This streamline strikes the disk at a 
radius sin^(6'o)-Rd = (1 — £^)Rd- 

Geometrically, a fraction e of the starlight strikes 
the cavity inner edge because the remaining 1 — e is 
cleared by the outflow. However, smaller e leads to 
a broader outflow cavity, which causes a greater por- 
tion of the reprocessed starlight to reach the disk. Af- 
ter performing ray-tracing calculations of reradiation 
from the inner edge of the infall, using a geometry 
like that shown in figure El we find that fp is ade- 
quately described by O.le"*^'^^ in the relevant range 
20% < e < 90%. The infall envelope is translucent to 
this reprocessed radiation, but our estimates indicate 
that shadowing of the disk by the infall is insufficient 
to change this result significantly. 

The disk's edge is too hot to fragment unless 



period > 2.4x10 





1.8 rp 


^10K\ 




4500 K 


^ Tc ) 



cos 00 



(42) 



9/4 

yr. 

(43) 

Consider now the evolution of the effect of irradia- 
tion as the accretion rate declines. The critical disk 
midplane temperature, eq. (|25|) . varies as mUc . The 
radiative flux normal to the outer disk varies as the 
central luminosity (oc Mace for low- mass stars); if the 
infall is optically thin then there is another factor 
of Mace in the reprocessed fraction. The flux there- 
fore varies as Ml^^c'i implying a midplane tempera- 
ture that varies as Maec^"*^'^- ^^^^ weaker dependence 
implies that a disk which is stabilized by irradiation 
will remain so as the mass supply wanes. This might 
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Fig. 2. — A schematic diagram for irradiation in the presence of an outflow cavity, in the theory of H4.2.2I The inflow envelope 
is excavated within a particular streamline (eq. 1421 . solid line); visible starlight is absorbed near this innermost streamline and 
reprocessed into infrared light that illuminates the disk (dashed arrows). 



not hold if accretion were to stop suddenly compared 
to the disk's viscous time, but this is not realistic. 

To summarize, local fragmentation due to Toomre 's 
instability is quenched by irradiation in protostellar 
disks, in all but exceptionally long-period systems 
whose orbital separations are of order 1000 AU. 

However, disks fragment at smaller periods if they 
form from cores that are warmer than 10 K or 
from more turbulent cores than were described in 
^ Disks around ma ssive stars fall in this category 
(|McKee k Tanll200,'^ . although their central illumi- 
nation is different than we have assumed. 

5. GLOBAL INSTABILITY 

In light of the previous section, one expects irradia- 
tion to prevent local fragmentation and possibly also 
to prevent the enhanced angular momentu m trans- 
port f rom self-gravity. Note, however, that iGammij 
(|2001l )'s results were for razor-thin disks, and proto- 
stellar disks are rather thick in comparison. There- 
fore, the disk mass can become large enough to induce 
global instability. As we shall see, the threshold for 
this depends on the ability of local angular momen- 
tum transport (e.g., due to MRI) to remove material 
fr om the disk. 
Adams et all (|l98fll ^ and lShu et iZI t99(i ) describe 



a global "SLING" instability in which the disk be- 
comes massive enough that the star is perturbed from 
the center of ma ss by the action of a lopsided spiral. 
Shu et aP (fl990") provide the following formula (their 
eq. [Ill]) for the criterion of instability: 



Md + M 47r 



MiQ) 



(44) 



where 



MiQ) ^ 



(45) 



where the approximation is exact in the limits Q — > 1 
and Q — > oo and holds within 2.3% for all Q > 1. 
When it sets in, this instability is thought to trans- 
port material rapidly eno ugh to prev ent fragmenta - 
tion ( Laughlin et al.l 19981 . see however Bonnelll994^ . 

We shall not posit Q = 1, because we have seen 
in ^4.21 that irradiation tends to warm the disk so 
that Q > 1. Instead, let us assume S oc r"^/^ 
so that Md = 47ri?^S. If we approximate 17^ ~ 
G{M + Md)/Rd and retain H = Cs/9., then ^ be- 
comes 

«-?(-^)- 

With equations (|44j) and (|45|) . this gives a maximum 
stable disk mass as a function of H/r. We find that 
Q > 1 and Md < M requires 0.06 < H/r < 0.28; 
within this range, we approximate the exact solution 
of equations (011), (HSJ, and 



^< 3.09- + 0.13 

M r 



(47) 



for global stability. The approximation is valid within 
1%. 

This can be compared to the mass accumulated, 
which for S oc r"^/^ is 



Md 



Mace 4 



r y 

IT 3a [hJ 



(48) 



The comparison sets a minimum for a/ag, where 



The theory of 321 indicates that 



-2.6 Jj 



J. 

3/4 



(49) 



(50) 
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independently of the physical scale of the parent 
cloud. 

Taken together, we find that the minimum a re- 
quired to keep the disk stable can be approximated 
(within 4.2%) 



a > 0.52 



H 



2.76 



(51) 



The instability is clearly very sensitive to the mid- 
plane temperature. Computing H/r = Cg/vk using 
the irradiation temperature model of ^4.2.21 (which 
gives H/r ~ 0.11 for fiducial parameters), we find 
that stability requires 



a > 0.19 



1.62 



fl.27 



60% \ 



0.27 



V lOK 



1.27 



165 Mq pc- 



0.4 



/4500\°''^6 

\~Tr) 



(52) 



It has been argued that local gravitational trans- 
port becomes significant when Q falls below ~ 1.4 
( Laughlin &: Bod enheimer 1994,). In its absence, the 
disk is left with magnetoturbulent (MRI) transport, 
which is thought to saturate at much lower values of 
a in a weakly ionized medium. 

With eq. (|^^ . this indicates that global spiral 
modes can be required for protostellar disks to process 
the mass supply, depending on the infall parameters 
and on the detailed irradiation temperature (which 
we have calculated only approximately). Disks that 
are especially thick or whose parent cores are espe- 
cially slowly rotating {fj « 1) will be stable at 
smaller values of a. 

6. SUMMARY AND CONCLUSIONS 

6.1. Summary 

We derived in ^ a simple set of relations for the 
smallest and most thermally supported collapsible 
objects, i.e., cores, that can exist at a given mean 
gas temperature in a parent cloud of a given column 
density. The cloud column sets both the scale for 
the cores' bounding pressure, and the coefficient in 
the line width-size relation for cloud turbulence. This 
turbulence aids in core support and increases the min- 
imum core mass, while also accelerating the accretion 
rate when cores do collapse f H2.2|) . In addition, tur- 
bulence bestows cores with rotational angular mo- 
mentum that determines the characteristic scales for 
protostellar accretion disks f H2.1() . This theory does 
not account for the role of a static magnetic field, 
nor does it prescribe the distribution of overpressures 
expected in a realistic model for molecular cloud tur- 
bulence; nevertheless, it matches the observations of 



prestellar cores of a couple solar masses f ^2.4() with a 
minimum of free parameters. In our model, substel- 
lar cores would derive from regions of especially low 
temperature and high degrees of overpressure. High 
mass cores could come from the opposite conditions 
or, more likely, derive the bulk of their initial support 
from nonthermal pressure. 

Turning next to the fragmentation boundary as 
predicted even in axisymmetric models of protostellar 
collapse, we found in ^that the cores are typically 
su bject to a phase of prompt frag mentation according 
to iMatsumoto Sz Hanawa criterion. This is 

a possible mechanism for the production of relatively 
tight binary stars. However, we note that this occurs 
before most of the core mass has accreted; therefore 
these fragments may be driven to merge by radiating 
angular momentum to the accret ion disk. We also 
stress that IMatsumo to Sz Hanawal 's simulations used 
an approximation to the thermal evolution of the col- 
lapsing core, and that we and others have found the 
gas temperature to be a controlling physical factor in 
fragmentation. We therefore do not consider it cer- 
tain that binaries form by this mechanism. 

Examining the local gravitational stability of proto- 
stellar disks in N4.ll we mad e extensive use of the sim- 
ulations bv iGammid (•2001') to identify a critical disk 
temperature that divides accreting from fragmenting 
disks. For disks heated solely and locally by their own 
viscosity, we found in H4. 1.11 that local fragmentation 
is impossible for disk periods shorter than about 460 
years (depending on the opacity law), and that frag- 
mentation is suppressed until much longer periods if 
the disk is optically thin to its own cooling radiation. 
This conclusi on is enti r ely co nsistent with the recent 
argument by iRafikovl ( 20041 ^ that the f ormation o f 
giant planets by direct disk instability ( Boss! 19971 1 
is impossible on thermodynamic grounds. Moreover, 
when we account for heating by the star's accretion 
luminosity (reprocessed by the infall envelope), we 
find that the local instability is entirely quenched out 
to very long periods, ~ 10^'^ years. A decline in the 
accretion rate does not alter this conclusion. 

Lastly, we consider the onset of global instability 
due to the accumulation of a fini te disk ma s s. Ap - 
plying the instability criterion of Shu et all l|l990l ). 
we find that the intrinsic viscosity parameter a (e.g., 
due to MRI) must exceed roughly 0.2 in order to pre- 
vent global instability. This makes global instability 
a primary candidate for angular momentum trans- 
port within protostellar disks. When it does occur, 
this instability is likely to saturate and provide rapid 
angula r momentum transpo rt rather than fragmen- 
tation ( Laughlin et al.lll99^ 1. 
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6.2. Implications for Numerical Simulations of 
Multiple Star Formation 

Thermal evolution plays a critical role in determin- 
ing th e onset and outcome of gravitational stability in 
disks flToomrelll964l:ICa,sseTi et a].lll98ll:lTnrnlev et al " 

" li 9991 . l2nnnl 



1^ 



2m 



h994l:IPickett et a]jll99Hl. 
IPickett et a1.ll2nn^. 



ammie 



We find that local disk 



fragmetation is quenched by stellar irradiation except 
at ve ry large periods (§ 14.21 see also IP'Alessi o et al 



1997t) and, even in the absence of irradiation, inhib- 
ited for short periods (< 460 years) by viscous heating 

These observations are quite important for the in- 
terpretation of large-scale numerical simulations de- 
signed to predict star formation. Currently, only a 
few multidimensional calculations have accounted fo r 

,i.l l2nnnh . 



radiation transfer during collapse (jB^ 
and then only in an approximate manner. No dynam- 
ical simulations have yet included disk irradiation due 
to starlight that is reprocessed by the infall envelope. 
Disk fragmentation in these simulations is unlikely to 
be realistic. 

For instance, three quarte rs of the b r own d warfs 

arise 



formed in the simulations of lBate et al. I (|2nn2l ' 



via disk fragmentation. Since these simulations em- 
ploy a barotropic equation of state, we predict that 
these brown dwarfs would not form if stellar irradia- 
tion were accounted for. 

We have demonstrated in ^that initial conditions 
for most collapsing cores are f avorable for the prompt 
fragmentation simulated by iMatsumoto &: Hanawal 
(|2003). However, since these were also performed 
with a barotropic equation of state, and since they 
do not extend long enough to determine whether the 
fragments will merge, we consider them provisional 
and await radiation transfer calculations of the same 
parameter space. 



6.3. Conclusions: Stellar Binaries and the Brown 
Dwarf Desert 

Our foremost conclusion is that stellar and substel- 
lar companions are not produced by disk fragmenta- 
tion, except possibly in an early phase of core col- 
lapse (although we argue this is unlikely: see or 
at very large periods of ~ 10^'^ years f ^4.2|) . This 
indicates that turbulent perturbations beyond rota- 
tion and compression are responsible for the pro- 
ducti on of binary stars ( Klein et al.ll2003l : iBate et al 
20o3). Since companions formed in disk instabili- 



ties would be initially of substellar mass, this closes 
a formation channel for substellar companions over 
a wide range of o rbital periods. Observations (e.g., 
McCarthv &: Zuckerman-2004 ') show a marked dearth 
of companions in the mass ran ge between stars and 
planet s - a brown dwarf desert. lArmitage &: Bonnelll 
(200i) have ascribed this to inward disk migration 
that destroys brown dwarfs. However this presup- 
poses that they can form within disks; here we have 

argued that they cannot.^ 

Barrado v Navascues e t al.l (|2004l ^ argue from their 
observations of disk accretion onto free-floating 
brown dwarfs that these must form in isolation rather 
than from fragmentation in a disk or in a collaps- 
ing turbulent core. Our theoretical results also argue 
against a disk origin for these objects. 
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from our anonymous referee have been very helpful in 
improving the clarity of the paper. CDM's research is 
funded by NSERC and the Canada Research Chairs 
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SPECIFIC ANGULAR MOMENTUM DISTRIBUTION OF PRESTELLAR CORES 

We now compute the angular momentum distribution quoted in equations (fTn|) and (fTT|) according to the 
model for core rotation described in ^2.11 Specifically, we assume in the present calculation a spatially homo- 
geneous velocity fields v whose components are uncorrelated Gaussian fields obeying a line width-size scaling 
like that in equation The effect of gas pressure, which introduces correlations by inhibiting compression, 
is discussed separately in ^2.11 As a consequence of homogeneity, we assume no correlation between v and the 
core density field p{r). See ^2.31 for a discussion of these and alternative assumptions. 

Since we are interested in the specific angular momentum normalized to Rat^T{R), we compute j and o"nt(-R) 
separately. 

Angular momentum 

To satisfy the line width-size relation in this model, we must impose that the velocity difference between the 
two points scales as a square root of the distance between the points. More precisely, 

{h{n)-v,{r2)]') = k\r,-r2\5,, (Al) 

where we assume that |ri — <C v'^/k if v = (Ivp)-'^/^ is the characteristic turbulent velocity. Here, angle 
brackets denote an ensemble average over realizations of the turbulent velocity field, or equivalently, over 
the location of the core within this field. Equation (jAll) is assumed to hold for each velocity component Vi 
independently. Expanding equation (|Aip . 



{vi{ri)vj{r2)) 



1 



1 



k\ri 



r2 



(A2) 



.3 2 

Skipping to the result, the root mean square angular momentum of a core in this velocity field is given by 

(l2)V2 = 27rk'/^R^/^po [F{p)f^ , (A3) 

where R is the core radius and po is a reference density; p(f) = p{r)/pQ is the spherically symmetric normalized 
density profile expressed in terms of the dimensionless radius f = r/R, and the form factor F is given by the 
following double integral: 



1 

^0 



dridr2p{ri)p{r2)rir2{r'l + r'lf'/'^ x 

l+.)5/2_n _ )5/2] 



1 r 
5 



(l + (?)3/2-(l-g) 



3/2' 



(A4) 



Here q = 2rir2/(rf + r^). 
Evaluated over a Bonnor-Ebert sphere, i^^/^ = 0.0294. Such spheres obey M = 0.735poR'^ , so that 



{j 



2\l/2 



(L2)1/2 



M 



0.251k^/'^R^^'^. 



The derivation is as follows. The total angular momentum is given by 



+ + L^. 



(A5) 



(A6) 



By spherical symmetry the components of angular momentum are uncorrelated; this can be checked directly. 
Therefore, 

{L') = 2>{Ll) = m,-L2f), (A7) 

where 



L2 



d'^rp{r)xvy, 
j d^rp{r)yvx 



(A8) 
(A9) 
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But {L1L2) = and, by spherical symmetry, (Lf) = {L^)- Therefore, 

(L2) = 6(L?) 



(AlO) 



6 



d^rid^V2 p{ri)p{r2)xiX2{vy{vi)vy{r2)) . 



We can now use Eq. ()A2|) and evaluate the integral above. First, note that only the second term on the 
right-hand side of Eq. HA2() contributes to the integral. Second, because of the spherical symmetry we can 
substitute X1X2 by (rir2)/3 in the integrand. Thus we have for the total angular momentum: 



(L^) = -k j j d?vid^r2p{ri)p{r2){YiY2)\vi - Y2\ 



(All) 



The integrand in Eq. HA11|) depends only on the magnitudes of ri and r2, and on the angle 9 between them. 
Therefore, the other three free variables are integrated out, and we are left with the 3-d integral: 

' 2 
2 



— k j dri4:Trrfp{ri) x 
J J dr 22Tid cos cos 6 p{r2)rir2\J r\ + r. 



\/l — q cos 9 



where 



2rir2 
r? -|- r. 



2 • 



(A12) 



(A13) 



The angular integration in Eq. HA12|) can be performed analytically, most simply by using ^/Y^^q'cx^9 as an 
integration variable. The result of this calculation is presented in Equations HA3|) and HA4|) . 

One- dimensional velocity dispersion 
We now compute the one-dimensional velocity dispersion: 



a 



j d^r ^[^^(r) - Vx] 
d^r p 



= j d^r pv^irf' 
where p = p/M and v denotes the center of mass velocity: 

V = I d^^Y pv{r) 



(A14) 



(A15) 



In equation (|A14I) . the first term reduces after ensemble averaging to f^/3, according to equation HA2|) . 
Likewise for the second term, 



d^rid^Y2 p{ti)p{t2){v.j:{yi)vx{t2)) 



(A16) 



d rid Y2 /5(ri)/5(r2) -v - -k\ri - r2 



1 



1 



.3 2 

The constant terms cancel in the ensemble average of eq. HA14|) . leaving 

1 



-k 



dpYiSr2 p(ri)p(r2)|ri - r2| 
^k j ATTrldrip{ri) J 27rr|p(r2)-^rf -|- x 
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J dcos9\/l — qcosO 



kR 



dridr2 rjr|/5(ri)p(r2)^rf + x 
(1 + qfl' - (1 - qf'^^ 



Evaluated over a Bonnor-Ebert sphere, 



and in combination with equation (jASI 



(a2)i/2 = o.607fei/2i?i/2^ 



(j2)V2 = 0.414(^72)1/2^. 



(A17) 



(A18) 



(A19) 



Since the velocity distribution is Gaussian, all of its linear moments (e.g., and a) are Gaussian distributed. 
The cumulative probability distribution of j, V{< j), is therefore Maxwellian with a width fixed by equation 

dV{<j) f ( 3j2 



d(jy(j2)i/2) 



2(f> 



(A20) 



A log-normal fit to this distribution is \ogiQ j / {p)^^"^ = —0.080 it 0.52, but this is not especially accurate. A 
more faithful representation is log^g ^/(j^)^''^ = — 0.088]'^Q4g. 



